Spontaneous Breakdown of Superhydrophobicity 
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In some cases water droplets can completely wet micro- structured super hydrophobic surfaces. The 
dynamics of this rapid process is analyzed by ultra-high-speed imaging. Depending on the scales of 
the micro-structure, the wetting fronts propagate smoothly and circularly or - more interestingly - 
in a stepwise manner, leading to a growing square-shaped wetted area: entering a new row perpen- 
dicular to the direction of front propagation takes milliseconds, whereas once this has happened, 
the row itself fills in microseconds ( "zipping"). Numerical simulations confirm this view and are in 
quantitative agreement with the experiments. 

PACS numbers: 83.50.Rp,68.03.Cd,05.20.Dd,02.70.Ns 



Micro-structured materials can show a superhydropho- 
bic behavior with effective contact angles of 160° and be- 
yond ("Lotus effect"). For many applications this effect 
is wanted: when a liquid droplet is deposited on micro- 
textured hydrophobic surfaces, it can bead off cornpletely 
and wash off contaminations very efficiently [11,0,1 
Superhydrophobic materials are now used for a wide se- 
ries of applications in medicine, coatings, self-cleaning, 
textiles, and microfluidics However, under cer- 

tain conditions, the superhydrophobicity ( "Cassie-Baxter 
state" [lo|, hereafter CB) spontaneously breaks down 
[ill , \v\ : fluid enters in between the micro-structures and 
spreads, resulting into a smaller contact angle ("Wen- 
zel state" [m, hereafter W). In many cases this tran- 
sition from the CB to the W state is highly desirable . 
For example, water-repellent layers on leaves keep plants 
healthy, but also complicate spraying crops effectively 
with pesticides [13]. Other examples are heterogeneous 
porous catalysts, where superhydrophobicity is an un- 
wanted effect as it reduces the contact area [15|. 

The objective of this Letter is to characterize and an- 
alyze the dynamics of the spontaneous breakdown of su- 
perhydrophobicity in terms of geometry, wetting, liquid- 
gas surface tension cjig^ and the Hquid viscosity rj. 

The fabrication of highly precise and controllable 
micro-structured surfaces becomes possible through a 
micro- molding technique [16]. As material, we choose a 
solution of styrene-butadiene-styrene (a block-copolymer 
commercially available as "kraton") dissolved in toluene. 
Kraton is used because thin layers are easily removable 
from the mold. The thickness of the resulting translucent 
film is only ~ 40 //m, allowing for optical imaging of the 
CB to W transition through the bottom of the film. 

For the experiments presented here we choose a reg- 
ular, periodic structure of square pillars with height 



h = 10/im, width w = (both kept fixed for all our 
experiments), and as control parameter a gap width a 
between 2 jam and 17 jim. The wavelength d = w ^ a oi 
the regular square lattice is thus between Ijim and 22/im 
(see figure [2]). 

Water droplets are deposited on the surface with a sy- 
ringe pump at a very low flow rate (5/il/min) from an 
outlet of a vertical thin tube (outer diameter 0.158 mm). 
The outlet is set parallel to the flat film, 2 mm above 
it. A typical droplet grows slowly and reaches the dry 
surface within a minute. Then the flow is stopped and 
a (meta-)stable drop in the CB state (contact area ~ 1 
mm^) is observed. For the chosen micro- textured sur- 
faces the (effective) contact angle in the CB state is ^ ~ 
160°, while for a smooth surface made of the same ma- 
terial it is l9 - 100°. 

At some point the CB state can spontaneously break 
down (figure [T]). 

The extremely rapid filling process towards the W state 
is recorded from the bottom of the film (see figure [3^ for 
a sketch) with a microscope and a high-speed imaging 
system (Photron Ultima APX-RS), using frame rates in 
the range between 4000-50000 fps. In figure [3] the large 
scale time evolution of the filling process is shown. After 
fluid sinks down at a certain point, a lateral spreading de- 
velops which sensitively depends on the parameters char- 
acterizing the micro-textured surface: while for large gap 
width a = 11 jim a round shape of the fully wetted area is 
emerging (figure [3h), at smaller a = 5/im the propagat- 
ing fronts reflect the structure of the underlying square 
lattice, leading to a square-shaped wetted area (figure 
[Sb). At even smaller a = 2/im we never observed any 
transition to the W state. 

For the wetting cases, the velocity of front propagation 
drastically depends on the gap width a (see figure |4]): For 
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FIG. 1: Transition from the (meta-) stable CB to the W state. 
A drop softly deposited on a micro-patterned surface can stay 
suspended with air pockets trapped in the grooves underneath 
the liquid (left). At some point the CB state spontaneously 
breaks down. The drop then homogeneously wets the sub- 
strate, resulting in a lower contact angle (right). 




FIG. 3: (color) Bottom views of the front evolution of the 
transition, as sketched in (a). In (b) three snapshots for the 
case with a = ^iim are shown, leading to square-shaped wet- 
ted area. In (c) it is a = ll/xm, resulting in a circular wetted 
area. Figures (d) and (e) show the results of the correspond- 
ing numerical simulations with the Lattice Boltzmann method 
with a — and ll/xm, respectively. The infiltration point 
has been centered in the figure but is not located in the center 
of the droplet 's base: its dependence on local energy barriers 
can indeed produce uncertainty in its location. 
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FIG. 2: Sketch (a) and scanning electron microscopy (SEM) 
picture (b) of the micro-patterned substrate. The geometrical 
dimensions of the pattern are the pillar height h — lO/xm, the 
pillar width w — 5fim (both kept fixed for all experiments), 
and the gap width a, which is varied between 2/jm and 17/im. 
The wavelength d of the pattern is d = a-\-w. In (b) we have 
a = ^iim. Typical equilibrium contact angles for a smooth 
surface made up of the same material are ^ 100°. 



a = 11 /im we measure a mean velocity v ~ 700mm/s, two 
orders of magnitude faster than the case with a = 5/im, 
where v ^ 7 mm/s. Moreover, in the latter case the 
front propagates in a stepwise manner (figure [4^) and 
fronts at different positions on the sample or in different 
directions show a remarkable variation in their velocity 
("dispersion"). 

For the more striking case of the square-shaped wetted 
area (a = 5/im) we show details of the wetting dynam- 
ics in figure [5^. The front is here advancing from top 
to bottom and is slowing down while wetting occurs in 
a zipping manner: the timescale of the front for entering 
a new row is slow as already shown above (figure UK), 
namely typically d/v ~ 1.4ms, the timescale for the fill- 
ing of the once entered row itself is about two orders of 
magnitude faster, ~ O.Olms. 

In order to better reveal the physics of the wetting 
mechanism, we performed numerical simulations with a 




FIG. 4: Wetting front evolution for micro-patterned surfaces 
with different gap-width a. We show a case where a = ^iim 
(a) and a case where a — 11 /jm (b). In both cases, nine 
different experiments were carried out with the same sample. 
The position of the front has been normalized with respect to 
the pattern wavelength d = a -\- w and time has been made 
dimensionless with to as in equation (|4|. Note the difference 
in front speeds v in the two cases: in the small wavelength 
case (left) a critical slowing down (v ^ 7 mm/s) induces large 
dispersion and zipping is revealed. The large wavelength case 
(right) is less dispersed and much faster, reaching velocities 
ol V ^ 7 • lO^mm/s. 



three-dimensional Lattice Boltzmann algorithm [12| for 
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FIG. 5: The lateral zipping mechanism (left to right) close 
to the critical point (0 ~ Oc). (a) Snapshots of the zip- 
ping mechanism recorded in an experiment with a = 5/jm 
with the front propagation from top to bottom, (b) The 
results are also compared with snapshots from numerical 
simulations, (c) Numerical results for the front propaga- 
tion (lines) in the lateral zipping process close to the criti- 
cal point. From top to bottom we show the following cases: 
cos 0/ cos Oc ^ 0.868, 0.890, 0.927, 0.945, 0.971. The geometri- 
cal aspect ratio for the simulations is chosen to be the same 
as in the experiment with a = 5/im, whose typical outcomes 
are also reported (o). 



single component multiphase fluids. Wetting proper- 
ties with surface tension are introduced as explained in 
[lit , leading to wetting angles comparable with those in 
the experiments. Geometrical structures are reproduced 
with the same aspect ratio as in the experiment. Sim- 
ilarly, in the large gap width case a = 11 fim (nearly) 
spherical wetted areas are observed (figure [3^) , whereas 
in the case a = 5/im we observe square-shaped wetted 
areas, see figure [3]i. Also the zipping wetting behavior is 
reproduced in these simulations for a = bfim (see figure 
When repeating the simulations with a gap width 
a = 2/im, we do not observe any lateral infiltration of the 
substrate and the CB state remains stable, again, just as 
experimentally observed. How to explain the transition 
from CB to W ? The energy of a droplet in equilibrium on 
a substrate is monotonically increasing with the effective 
contact angle [lH. As a result, when the CB effective 
contact angle is higher than the W one, we would ar- 
gue that the CB droplet always collapses towards the W 
state. Anyhow, intermediate states with higher energy 
can be encountered [lH: they represent energy barriers 
to overcome and work has to be done in order to induce 
the transition 0, [13] • The bottom- view observations re- 
veal that the micro-structured surface is not wetted all 



at once, but the infiltration starts locally. In our case, 
this happens spontaneously after a few seconds whereas 
in other cases (i.e. larger h) the infiltration can take 
longer and one may want to trigger it. Once the tran- 
sition has started, how to theoretically understand the 
shape of the wetted area, the zipping, and the different 
involved timescales ? For the dynamics of liquid moving 
in between the posts, both the free surface and wall in- 
teractions play a role. A small advancing displacement 
<C a of the interface within the posts would take place 
in favour of a gain in surface energy due to the reduction 
of liquid-gas interface {aigadx) on the top of the pillars. 
This is a pulling mechanism for the transition that must 
be balanced with the energy loss to wet a small portion 
of hydrophobic wall {aig {2h + a)dx cos 6>, cos ^ < 0). The 
overall energy gain is 



dEs = aigadx + crig{2h + a)dx cos 0. 



(1) 



The limit dEg characterizes a critical contact angle 
Oc for the wetting properties above which the liquid can 
spread horizontally through the posts. For square posts 
and square geometry it is equal to 
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Note that in this equation there is no reference to the 
width w of the posts because we describe the horizontal 
filling, different from a simultaneous vertical collapse [l|, 
[l2[. With this critical value, the energy gain ^ can be 
rewritten as 



dEs = (Tigadx I 1 



cos 
cos 6 c 



(3) 



where for cosO < cos Oc no favorable propagation is ex- 
pected while in the other case, when cos 6c < cos 6 < 0, 
the propagation is energetically favored. For the two an- 
alyzed cases a = 5/im and 11/im we obtain Oc = 101.54° 
and 110.78°, respectively. As ^ = 100° for water on flat 
kraton, we thus understand that for gap widths smaller 
than a = 4. 2/im the propagation to the W state is ener- 
getically not possible. 

We can also understand the critical slowing down of 
the front close to the critical angle Oc (figure [4j). There- 
fore we have to identify the relevant time scale r for the 
system and quantify its fluctuations with respect to the 
geometry. To do so, we estimate the energy cost in terms 
of dissipation and compare it with the change in surface 
energy of equation (|3]). With the small dimensions un- 
der consideration, flows can be regarded as laminar and 
viscous dissipation is dominating as compared to flow 
inertia. A simple estimate for the rate of viscous dissipa- 
tion per unit volume in a fluid with viscosity r] is e ^ 7^7^ 
[Tot , with 7 ~ 1/r the characteristic shear rate. If we 
integrate it over the small volume dV = ahdx and time 
lag r, we obtain the viscous contribution for the energy 
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as dEd ^ erdV = rjahdx/r. Balancing dEg in equation 
(|3]) with dEd, we deduce the time scale 



to 



(1 



cos 9 

cos 



(4) 



with to = Tjh/cFig ^ 10~^5 for typical values of 7^, a^^, and 

= 10/im. Equation (|3|) reveals the origin of the criti- 
cal slowing down of the front propagation, its dispersion 
and the appearance of zipping: if ^ ~ 6>c, the timescale r 
of front propagation diverges as r ~ l/\0 — Oc\' Around 
this critical pointy a small uncontrollable variation (im- 
perfections of the microstructure, dust deposit) in the 
local wetting angle translates into a huge dispersion in 
the time scales, which is consistent with the experimen- 
tal observation shown in figure 2^. A flat front is slowly 
allowed to proceed further, and entering a new row per- 
pendicular to the direction of the front can take up to 
milliseconds, whereas once this has happened, the row it- 
self is filled on a faster time scale {zipping). In this way a 
square-shaped propagating wetting pattern emerges (see 
the case with a = 5/im in figure [3)3 where Oc ^ 101.54° 
and - 100°). 

On the contrary, when approaching the other limit 
cosO/ cos Oc <C 1 (as in figure [3b, where ~ 100° and 
a = 11/im, implying Oc ^ 110.78°) the dynamics becomes 
more and more determined by the time scale to itself (see 
figure [4)3) and the propagation through the posts is ex- 
pected to be smooth. In this limit it is the pure pulling 
mechanism of surface tension in equation ([3]) that domi- 
nates the spreading dynamics and is prevalent over wall 
effects [20|. For this reason, the geometrical properties 
of the lattice do not emerge and the front assumes an 
almost circular shape. 

In conclusion, we have experimentally, numerically, 
and theoretically revealed the origin of zipping wet- 
ting behavior at the spontaneous breakdown of super- 
hydrophobicity. We observed that the wetting process 
starts locally from a single point and then proceeds lat- 
erally, depending on a critical contact angle. 

Close to this critical point the driving energy for the 
lateral filling reduces, the front propagation slows down 
through viscous effects and zipping wetting is observed. 
As a consequence, the front propagates in a stepwise 
manner and square-shaped [2l| wetted areas emerge. Al- 
though we have observed the origin of this zipping pro- 
cess, its dynamical details have not been characterized 
here and will be the object of a forthcoming paper [22]. 

The critical contact angle (equation ([2])) can be cal- 
culated from the geometrical properties of the micro- 
textured material and is consistent with our experimental 
and numerical findings. 

Our results are useful for the design and fabrication 
of micro-structured surfaces with certain wetting proper- 
ties: for given hydrophobic material and liquid, the post 
height and the gap width can simply be calculated to 



obtain some desired critical contact angle. 

The next step will be to control the time scale of the lo- 
cal initiation point, characterize the dynamics of the zip- 
ping process and extend our analysis to more complicated 
micro-patterned structures, including three-dimensional 
ones . 
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